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1 Introduction 

Events with high transverse momentum jets in proton-proton collisions are described by quan- 
tum chromodynamics (QCD) in terms of parton-parton scattering, where the outgoing scat- 
tered partons manifest themselves as hadronic jets. Measurements of the inclusive jet and dijet 
cross sections can be used to test the predictions of perturbative QCD, constrain parton dis- 
tribution functions (PDFs) of the proton, differentiate among PDF sets, and look for possible 
deviations from the standard model. 

In this paper, measurements of the double-differential inclusive jet (p + p — > jet + X) and dijet 
(p + p — > jet + jet + X) production cross sections are reported as functions of jet rapidity y and 
either jet transverse momentum pj or dijet invariant mass M«, at \/s = 7TeV. The data were 
collected with the Compact Muon Solenoid (CMS) detector at the CERN Large Hadron Collider 
(LHC) during the 2011 run and correspond to an integrated luminosity of 5.0 fb _1 , two orders of 
magnitude larger than the published LHC results from the 2010 run HHS- Jets are reconstructed 
up to rapidity 2.5, transverse momentum 2 TeV, and dijet invariant mass 5 TeV. The measured 
cross sections are corrected for detector effects and compared to the next-to-leading-order QCD 
predictions. 

2 Apparatus 

The CMS coordinate system has its origin at the center of the detector, with the z axis pointing 
along the direction of the counterclockwise beam. The azimuthal angle is denoted as (p, the 
polar angle as 9, and the pseudorapidity is defined as rj = — In [tan (6/2)]. The central feature 
of the CMS detector is a superconducting solenoid, of 6 m internal diameter, that produces an 
axial magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a lead 
tungstate crystal electromagnetic calorimeter (ECAL) and a brass/ plastic scintillator hadronic 
calorimeter. Outside the field volume and in the forward region (3 < \t]\ < 5) is an iron/ quartz- 
fiber hadronic calorimeter. Muons are measured in gas ionization detectors embedded in the 
steel return yoke outside the solenoid, in the pseudorapidity range \rj\ < 2.4. A detailed de- 
scription of the CMS apparatus can be found in Ref. QJ. 

3 Jet reconstruction 

The rapidity y and the transverse momentum pj of a jet with energy E and momentum p = 
(p x ,Py,Pz) are defined as y = (1/2) In [(E + p z )/ (E — p z )\ and p T = Jpl + pj. The inputs to 
the jet clustering algorithm are the four-momentum vectors of the reconstructed particle can- 
didates. Each such candidate is constructed using the particle-flow technique |5|, which com- 
bines the information from several subdetectors and is calibrated to account for the nonlinear 
and nonuniform response of the CMS calorimetric system to hadrons. Jets are reconstructed 
using the anti-fcx clustering algorithm [6j with distance parameter R = 0.7. The clustering is 
performed using four-momentum summation with the FASTJET package (7|, where the cho- 
sen distance parameter allows for the capture of most of the parton shower and improves 
the dijet mass resolution with respect to smaller sizes. The total transverse energy E Ej and 
missing transverse energy E™ 1SS are used in the event selection and are derived from the re- 
constructed particle-flow objects. They are defined as E Ej = Ejj, with Ejj = E, sin 0,, and 
gmiss _ _ v\ (£. s j n q { cos <^.£ _|_ £ { s j n q { s j n <£,y ), where the sum refers to all particle candidates 
and x, y are the unit vectors in the direction of the x and y axes. 
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4 Data samples and event selection 



The reconstructed jets require a small additional energy correction, mostly due to thresholds 
on reconstructed tracks and clusters in the particle-flow algorithm, and various reconstruc- 
tion inefficiencies. These jet energy corrections are derived using simulated events, generated 
by PYTHIA6 (version 6.4.22) |8l and processed through the CMS detector simulation based on 
Geant4 |9[, and in situ measurements with dijet, photon+jet, and Z+jet events Bui . These 
jet energy corrections correct reconstructed jets to the hadron level, as opposed to the par- 
ton level. An offset correction is also applied to account for the extra energy from additional 
proton-proton interactions within the same or neighboring bunch crossings (in-time and out- 
of-time pileup) ||T0| . The pileup effects are important for the lowest-px jets (10% jet energy 
scale correction and 1% systematic uncertainty [TlT| for jets with pj ~ 100 GeV) and progres- 
sively decrease with jet pj. For jets with pj > 200 GeV the pileup effects are negligible. The 
jet energy correction depends on rj and pj of the jet, and is applied as a multiplicative factor to 
the jet four-momentum vector. The factor is typically between 1.0 and 1.2 and is approximately 
uniform in rj. For a jet pj = 100 GeV the factor is 1.1, decreasing towards 1.0 with increasing 
pj. The typical jet pj resolution is 10% at pj = 100 GeV. The dijet mass Mjj is calculated from 
the corrected four-momentum vectors of the two jets with the highest pj (leading jets). The 
relative dijet-mass resolution, estimated from the simulation, ranges from 7% at Mjj = 0.2 TeV 
to 3% at Mjj = 3TeV. 

4 Data samples and event selection 

The data samples used for this measurement were collected with single-jet high-level triggers 
(HLT) H2| that require at least one jet in the event to have p T > 60, 110, 190, 240, or 370 GeV, 
respectively, in corrected jet transverse momentum. The online jet reconstruction uses only 
calorimetric information and the resulting HLT jets typically have worse energy resolution than 
the offline particle-flow jets. The lower-pj triggers were prescaled and the corresponding inte- 
grated luminosity of each trigger sample, £ e ff, is listed in Table [l] In the offline analysis, events 
are required to have at least one well reconstructed proton-proton interaction vertex |fT3| . In 
order to suppress nonphysical jets, i.e., jets resulting from noise in the electromagnetic and/ or 
hadronic calorimeters, the jets are required to satisfy the following identification criteria. Each 
jet should contain at least two particles, one of which is a charged hadron, and the jet energy 
fraction carried by neutral hadrons and photons should be less than 90%. These criteria have 
an efficiency greater than 99% for physical jets. 



Table 1: The integrated luminosity for each of the data samples. 



Min. jet trigger pj 
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190 


240 


370 


(GeV) 












Jet trigger name 


Jet60 


JetllO 


Jetl90 


Jet240 


Jet370 


^(pb- 1 ) 


0.41 


7.3 


152 


512 


4980 



The inclusive single-jet cross section measurements are made in five rapidity regions of size 
A\y\ = 0.5 over the range 0.0-2.5. Jets that do not satisfy the jet identification criteria are 
discarded and events are required to contain at least one jet that satisfies these criteria. In order 
to avoid any trigger bias, the jets are additionally required to have pj > 110, 200, 300, 360, and 
510 GeV for the five single-jet HLT triggers used, respectively. Figure [I] (left) shows the trigger 
efficiency as a function of the jet pj, for the central rapidity bin |y| < 0.5 and for the highest 
trigger threshold. The efficiency of each trigger path has been measured using events collected 
with a lower threshold single-jet trigger and confirmed with events collected with single-muon 
triggers. 
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Jets from rare events with instrumental noise, beam halo effects, and spurious signals that 
might survive the jet identification criteria are suppressed by requiring E^ lss /EEx < 0.3. Hard 
QCD processes do not generate true £™ 1SS and so the measured values of E™ 1SS in such events, 
due to jet mismeasurements and the energy resolution of the detector, are small compared to 
the total transverse energy. Hence, the distribution of the variable E™ 1SS /Y,Ej peaks close to 
zero for QCD events, while events contaminated with noise give larger values. Figure [2] (left) 
shows a typical distribution of the variable E™ lss /EEx for events with at least one jet with 
pj > 510 GeV, collected with the 370 GeV single-jet trigger. The data points exceed the QCD 
predictions at large values of E™ 1SS due to instrumental backgrounds and physics processes 
such as Z+jet(s), where the Z boson decays to neutrinos, and W+jet(s), where the W boson 
decays to leptons. 

For the dijet measurement, at least two jets with pj 1 > 60 GeV and pj 2 > 30 GeV and satisfying 
the tight identification criteria are required. If either of the leading jets fails the identification 
criteria, the event is discarded. The dijet measurement is performed in five rapidity regions, 
defined by the maximum absolute rapidity | y | max = max 1 3/2 1 ) of the two leading jets 
in the event. The use of the variable |y| m ax divides the phase space of the dijet system into 
exclusive rapidity bins, which correspond to different scattering angles in the center-of-mass 
frame. 

Low values of |y| m ax probe s-channel scattering at large angles, while large values of 1 1/ 1 max 
probe i-channel scattering at small angles. For each rapidity bin the trigger efficiency is ex- 
pressed as a function of My, and the events are required to satisfy a minimum mass threshold, 
which increases with |y| m ax- Figure [I] (right) shows the trigger efficiency for the central rapidity 
bin and for the highest trigger threshold. Because of the dijet topology, no further cut is needed 
on the E T 11SS /E £7 variable, as can be seen in Fig. (right). 
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Figure 1: Trigger efficiency as a function of the jet pj (left) and dijet mass Mjj (right) for the 
370 GeV single-jet trigger and for the central rapidity bins. 



5 Measurement of the differential jet and dijet cross sections 

In this section the construction of the jet transverse momentum and dijet mass spectra from 
the different samples is presented. Then the unfolding procedure, which translates the recon- 
structed spectra into true spectra, is described. Finally, the experimental uncertainties related 
to the measurements are described and discussed. 

5.1 Determination of transverse momentum and dijet mass spectra 

The jet pj (dijet invariant mass) spectrum is obtained by populating each bin with the number 
of jets (events) collected using the highest threshold trigger which gave more than 99% trigger 
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5 Measurement of the differential jet and dijet cross sections 
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Figure 2: Distribution of E™ 1SS /E £j for data events (black points) and simulated QCD events 
(continuous line) with at least one jet with pj > 510 GeV (left) and for dijet events with |y | m ax < 
2.5 and Mjj > 910 GeV (right), collected with the 370 GeV single-jet trigger. The distribution 
from the inclusive-jet selection is shown before the offline selection £^ 1SS /E£x < 0.3. The 
larger tail in the data is caused by other physics processes with true E™ 1SS (such as Z+jet(s), 
where the Z boson decays to neutrinos, and W+jet(s), where the Wboson decays to leptons), 
and instrumental backgrounds. 
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Figure 3: Spectrum construction from individual trigger paths. Left: inclusive jet pj spectrum 
for |y| < 0.5. Right: dijet mass spectrum for |y| max < 0.5. The different markers indicate 
different trigger paths. 



efficiency. Then, the yields from each trigger path are scaled according to the corresponding 
prescale value for this path (effective luminosity), as shown in Table [l] Figure [3] shows the 
reconstructed spectra for inclusive jets (left) and dijets (right), in the central rapidity bin, de- 
composed into the five contributing trigger paths. 

The observed inclusive jet yields are transformed into double-differential cross sections as fol- 
lows: 

dV = 1 Njets 

dp T d|y| e-£ eff Ap T A|y|' 

where Nj e t s is the number of jets in the bin, £ e g is the integrated luminosity of the data sample 
from which the events are taken, e is the product of the trigger and event selection efficiencies, 
both of which are greater than 99%, and Apj and A\y\ are the transverse momentum and ra- 
pidity bin widths, respectively. The width of the pj bins is proportional to the pj resolution 
and so increases with pj. The statistical uncertainty assigned to each pj bin takes into account 
the number of independent events that contribute least one jet in the bin. The largest fraction 
(more than 90%) of the observed jets in each pj bin originate from different events, however 
a small fraction of events contribute more than one jet. Such events are typically back-to-back 
dijet events, so closely balanced in pj that both jets end up in the same |y| and pj bin. 



5.2 Unfolding 
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The statistical uncertainty in the number of jets in a bin is e sta t = y/ (4 — 3/)/ (2 — /) • WNj e ts, 
where f = N\/ N ev is the fraction of events that contribute one jet in the given bin. The formula 
is valid under the assumption that the number of events that contribute more than two jets in 
each bin is negligible, which has been verified for the current measurement. 



The observed dijet yields are transformed into double-differential cross sections as follows: 



dV 



N 



dMjdlyl, 



r 



(2) 



where AMjj and A 



(ax are the mass and rapidity bin widths, respectively. The size of the dijet 
mass bins is approximately equal to or larger than the mass resolution at the bin center, while 
the bins at the edge of the spectrum have been merged to assure a minimal number of events 
in each bin. 



5.2 Unfolding 

Because of the detector resolution and the steeply falling spectra, the measured differential 
cross sections are smeared with respect to the particle-level cross sections. Each pj and mass 
bin contains events that have migrated in from neighboring bins and is missing events that 
have migrated out. For a steeply falling spectrum more events migrate into a bin than out. 
In order to allow for a direct comparison of experimental measurements with corresponding 
results from other experiments and with QCD predictions, the spectra are unfolded in order 
to correct for detector effects. The response matrix is obtained from the detector simulation 
and corrected for the measured differences in the resolution between data and simulation ||T0|. 
Figure |4]shows the response matrices for the jet pj (left) and the dijet mass (right) in the central 
rapidity bins. The unfolding is done with the RooUnfold package [14] using the D'Agostini 
method IfTHl . 
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Figure 4: Response matrices for the inclusive jet pj spectrum (left) and the dijet mass spectrum 
(right) in the central rapidity bins. 



5.3 Experimental uncertainties 

The dominant experimental uncertainties are related to the jet energy scale (JES), the luminos- 
ity, and the jet pj resolution. Other sources of systematic uncertainty, such as the jet angular 
resolution, are negligible. The agreement of the results for positive and negative rapidities has 
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5 Measurement of the differential jet and dijet cross sections 



also been confirmed. Figure [5]shows the effects of the experimental uncertainties in all rapidity 
bins for the cross section measurements. For rapidities up to \y\ = 1.5 the total uncertainty of 
both cross sections ranges from 5% at low pj or Mjj to 20% at high pj or Mp respectively. For 
higher rapidities the total uncertainty increases to 10-30% in both cases, with the exception of 
the highest dijet mass bin in the outer rapidity region of 2.0 < \y\max < 2.5, where the uncer- 
tainty is substantially larger. A discussion of the individual contributions to the uncertainty 
follows. 

5.3.1 Jet energy scale (JES) uncertainty 

The jet energy scale is the dominant source of systematic uncertainty. Because of the steep slope 
of the pj spectrum, a small uncertainty in the pj scale translates into a large uncertainty in the 
cross section for a given value of pj. The jet energy scale uncertainty is dependent on pj and tj 
and has been estimated to be 2.0-2.5% |rT|. The individual, uncorrected contributions to the 
JES uncertainty have been estimated and are discussed below. 

The JES uncertainty sources account for the pj and t] dependence of the JES within the total 
uncertainty. For the phase space of jets considered here, 16 mutually uncorrelated sources con- 
tribute to the total uncertainty, where each such source represents a signed lcr variation from 
a given systematic effect for each point in (pj, tj). Summing up separately the positive or neg- 
ative variations of the sources in quadrature will reproduce the total upward and downward 
JES uncertainties at each point. The uncertainties from all 16 independent sources are included 
in the HEPDATA record for this paper. 

The uncertainty sources are divided into four broad categories: pileup effects, relative calibra- 
tion of jet energy scale versus 77, absolute energy scale including pj dependence, and differences 
in quark- and gluon-initiated jets. The first category, containing pileup effects, has relatively 
little impact on the analyses presented in this paper. 

The second category, containing f/-dependent effects, parameterizes the possible relative varia- 
tions in JES, which for the dijet and inclusive jet analyses lead to correlations between rapidity 
bins. In principle these effects could also have a pj dependence, but systematic studies on 
data and Monte Carlo (MC) events indicate that the p T and t] dependence of the uncertainties 
factorize to a good approximation. 

The third category deals with the uncertainty in the absolute energy scale and its pj depen- 
dence and is the most relevant one for these analyses. The photon+jet and Z+jet events only 
constrain the JES directly in a limited jet pj range of about 30-600 GeV, and the response at 
higher (and lower) pj is estimated by MC simulation. The ^-dependent uncertainty arising 
from modeling of the underlying event and jet fragmentation is obtained by comparing pre- 
dictions from PYTHIA6 and HERWIG++. Most studies show that both generators agree with 
the data with differences comparable to those seen between data and MC. The uncertainty 
arising from the calorimeter response to single hadrons is estimated by varying the response 
parameterization by ±3% as compared to the central value. The final uncertainty arises from 
differences in the JES for quark- and gluon-initiated jets and is determined from MC studies. 

5.3.2 Luminosity uncertainty 

The luminosity uncertainty is estimated to be 2.2% |T6|. which can be directly translated into a 
2.2% uncertainty on the cross section normalization. It is fully correlated across all pj and mass 
bins. 



5.3 Experimental uncertainties 
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Figure 5: Effect of the relative experimental uncertainties for the inclusive jet (left column) and 
dijet (right column) cross section measurements, and for all five \y\ and |y| m ax bins, respectively. 
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6 Theoretical predictions 



5.3.3 Unfolding uncertainty 

The unfolding correction is closely related to the dependence on the pj and Mjj resolution and 
the spectrum slope. For the inclusive jet pj spectrum it varies between 5% and 10%, while for 
the dijet mass spectrum it ranges between 2% and 5%. The shape of the unfolding correction 
and uncertainty as displayed in Fig. [5] is understood as follows: the resolution in the observ- 
able, pj or My, improves when going from low to high values. As a consequence the effect of 
smearing is more pronounced in the lower pj or Mjj region. On the other hand the pj and Mjj 
spectra become steeper when approaching the kinematic limit at high pj or Mjj, leading again 
to a larger smearing effect than observed at medium values. 

The uncertainty introduced by the unfolding is caused by the modeling of the jet pj (dijet mass) 
resolution and the jet pj (dijet mass) spectrum in the simulation. In order to estimate the sensi- 
tivity of the correction to these inputs, the jet pj resolution is varied by ±10% and the jet (dijet 
mass) spectrum slope by ±5%. The former is motivated by the observed difference between 
data and simulation in the jet energy resolution [10J, and the latter is a conservative estimate 
based on comparisons of the theoretical and measured spectrum shapes. An additional con- 
stant 2% uncertainty is assigned to the dependence on the unfolding method. Overall, the 
unfolding uncertainty is of the order of 3^4%, fully correlated across the pj and mass bins. 

5.3.4 Other uncertainty sources 

The contributions from small trigger and jet identification inefficiencies, time dependence of 
the jet pj resolution, and uncertainty on the trigger prescale factor have been shown to be 
much smaller than 1%. To account for these residual effects a conservative uncertainty of 1% is 
assigned to each jet pj and dijet mass bin, uncorrected across the bins. 

6 Theoretical predictions 

The theoretical predictions for the jet cross sections consist of a next-to-leading-order (NLO) 
QCD calculation and a nonperturbative (NP) correction to account for the multiparton interac- 
tions (MPI) and hadronization effects. 

6.1 NLO calculations 

The NLO calculations are performed using the NLOJet++ program (v2.0.1) fT7\ within the 
framework of the fastNLO package (vl.4) ||18|. The renormalization and factorization scales 
(^r and }If ) for the inclusive and dijet measurements are identified with the jet pj and the 
average transverse momentum py ve of the two jets, respectively. The NLO calculation is per- 
formed using five different PDF sets: CT10 |H, MSTW2008NLO HO, NNPDF2.1 EU, HER- 
APDF1.5 [22J, and ABKM09 [23 1 at the corresponding default values of the strong coupling 
constant a s (M z ) = 0.1180,0.120, 0.119, 0.1176, and 0.1179, respectively. 

6.2 Systematic uncertainties 

The PDF variation introduces uncertainties in the theoretical prediction of up to 30%, while 
the variation of ocs(Mz) by ±0.001 introduces an additional 1-2% uncertainty. The uncertainty 
due to the choice of factorization and renormalization scales is estimated as the maximum 
deviation at the six points (ji F I >,?W » = (0.5,0.5), (2,2), (1,0.5), (1,2), (0.5,1), (2,1), where 
y. = pj (inclusive) or }i = p T ve (dijet). An additional uncertainty of at most 10% is caused by 
the nonperturbative correction. The scale uncertainty ranges from 5% to 10% for \y\ < 1.5 but 
increases to 40% for the outer \y \ bins and for high dijet masses and jet pj. Overall, the PDF 



6.3 Nonperturbative corrections 
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uncertainty is dominant for the high pj and high dijet mass regions. Figure |6]shows the effect of 
the systematic uncertainties for the two observables in all rapidity bins and for the NNPDF2.1 
PDF set. 

6.3 Nonperturbative corrections 

The nonperturbative effects are estimated from the simulation, using the event generators 
PYTHIA6 (tune Z2) and HERWIG++ 2.4.2 fl24|. (The PYTHIA6 Z2 tune is identical to the Zl tune 
described in [25J except that Z2 uses the CTEQ6L PDF while Zl uses CTEQ5L.) These models 
are representative of the possible values of the nonperturbative corrections, due to their dif- 
ferent physics descriptions. The nonperturbative correction is defined as the ratio of the cross 
section predicted with the nominal generator settings divided by the cross section predicted 
with the MPI and hadronization switched off. The central value of the nonperturbative correc- 
tion is calculated from the average of the two models considered, and ranges from 1% to 20%, 
being larger in the dijet spectrum because of the involvement of lower pj jets. 

7 Results 

The unfolded inclusive jet and dijet spectra are shown in Fig. |7j compared to the theoretical 
predictions. 

To compare the CMS data and the theoretical prediction, the ratio of the two is taken. Fig- 
ures |H]-[9] show this ratio using the central value of the NNPDF2.1 PDF set, accompanied by the 
total experimental and theoretical uncertainties. The theoretical uncertainties vary consider- 
ably among the different PDF sets, and in particular in the high-px arid high-Mjj region. The 
experimental uncertainty is comparable to the theoretical uncertainty. The additional curves 
represent the ratio of the central values of the other PDF sets to NNPDF2.1. Agreement is 
observed between data and theory in all rapidity bins, given the statistical and systematic un- 
certainties, with the various theoretical predictions showing differences of the order of 10%. 

8 Summary 

Measurements of the double-differential inclusive jet and dijet cross sections are presented us- 
ing 5.0 fb -1 of data collected with the CMS detector in proton-proton collisions at s/s = 7TeV. 
The measurements cover the jet pj range from 0.1 TeV to 2TeV, and the dijet-mass range from 
0.3 TeV to 5 TeV in five rapidity bins up to |y| = 2.5. The measured cross sections agree with 
the predictions of perturbative QCD at next-to-leading order obtained with five different PDF 
sets. Theoretical and experimental uncertainties are comparable, even at the limits of the ex- 
perimental phase space, so these results may be used to constrain global PDF fits. 
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Figure 8: Ratio of inclusive jet (left) and dijet (right) cross sections to the theoretical prediction 
using the central value of the NNPDF2.1 PDF set for the first three \y\ and |i/|max bins respec- 
tively. The solid histograms show the ratio of the cross sections calculated with the other PDF 
sets to that calculated with NNPDF2.1. The experimental and theoretical systematic uncertain- 
ties are represented by the continuous and hatched bands, respectively. 
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Figure 9: Ratio of inclusive jet (left) and dijet (right) cross sections to the theoretical prediction 
using the central value of the NNPDF2. 1 PDF set for the last two | y \ and | y \ max bins respectively. 
The solid histograms show the ratio of the cross sections calculated with the other PDF sets to 
that calculated with NNPDF2.1. The experimental and theoretical systematic uncertainties are 
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